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1. Introduction. 

Yang Mills theory on an l d continuum torus (d > 2) exhibits a phenomenon referred to as 
continuum reduction whereby the theory for I > l\ > is independent of I 1 fll], §]. At 
I = h, the theory goes from the confined phase (Oc: I > Zi) to the deconfined phase (lc: 
I < h). The order parameter is the Polyakov loop and rotational symmetry is spontaneously 
broken. More phases were conjectured to exist || in the continuum theory and these are 
refereed to as Xc phase with X = 2, • • ■ , d. 

The aim of this paper is to numerically establish the existence of the 2c and 3c phase 
in addition to the Oc and lc phase for the continuum Yang-Mills theory on a periodic torus. 
We will use the Polyakov loop to define an order parameter to be labeled, P [||, and it 
will take values in the range [0,0.5]. If the U(l) symmetry 2 under which the Polyakov loop 
transforms non-trivially is spontaneously broken, then P < 0.5. Let P x ,y z, be the order 
parameters in the three directions. Then, P x = P y = P z = 0.5 in the Oc phase. 

There are three possibilities for the lc phase and one of them is characterized by 
Py = Pz = 0.5 and P x < 0.5 with rotational symmetry still present in the (y, z) plane. It 
is difficult to numerically establish the order of the transition from the 0c to lc phase and 
we will leave it unresolved in this paper. 

x We use l\ to denote the critical size as opposed to l c since we will have a sequence of critical sizes. 
2 The U(l) symmetry is the limit of the Zn symmetry as N — > oo. 



- 1 - 



The 2c phase also has three possibilities and one of them is characterized by P x = P y < 
0.5 and P z = 0.5 with rotational symmetry present in the (x, y) plane. This transition 
occurs when I = fa < l\. One can argue that the lc to 2c phase transition is first order as 
follows. In the lc phase, P x was less than 0.5 and P y was 0.5. Since P x = P y in the 2c phase, 
it is necessary for at least P x or P y to change discontinuously at the lc to 2c transition. If 
one operator shows a discontinuity, all operators will generically show discontinuities and 
this signals a first order transition. 

Rotational symmetry is restored in the 3c phase and P x = P y = P z < 0.5. For the 
same reason as above, we expect the 2c to 3c phase transition to be of first order. This 
transition occurs when I = fa < fa- 

The 2c phase is characterized by two short directions and one infinitely long direction 
since the theory will not depend on the length of the direction where the U(l) is not broken. 
Therefore, this phase describes large N QCD in a small box at zero temperature (or infinite 
time). Confinement cannot be addressed in the 2c phase since we do not have large Wilson 
loops. The 3c phase describes large N QCD in a small box at high temperatures. The 2c 
to 3c transition is like the transition seen in perturbation theory on S 2 x S^Q where S 2 
replaces the two torus along which the U(l) symmetry is broken in the 2c phase. 

We extend our discussion to include l x x l y x fa torus with l x < l y < fa. The transition 
from 0c to lc will occur at l x = l\ and this is independent of l y and l z . The transition 
from lc to 2c will occur at l y = fa(l x ) with < l x < fa. Furthermore, fa(fa) = fa and our 
numerical results will show that £2(0) > 0. Finally, there is no dependence on l z . Continuing 
along the same lines, we can say that the 2c to 3c transition occurs at l z = fa(l x ,l y ) with 
^3 = fa(fajfa)- It is possible that one can obtain this critical size for small l x and l y by 
perturbation theory but it is necessary to consider the zero momentum modes of the gauge 
fields on all three directions. We do not address this problem in the paper. 

The results in this paper complement the results in the closely related paper by Bursa 
and Teper @j. We mainly focus on the continuum limit of the various phases using Polyakov 
loops. The paper is organized as follows. We define the relevant technical details in 
section ||. Our numerical results showing the existence of the various phases are presented 
in section |||. 

2. Technical details 
2.1 Lattice gauge action 

We used the single plaquette Wilson action given by 



n is a three component integer vector labeling the site, i labels a direction and i denotes a 
unit vector in the i direction. The link matrices U^n) are in SU(N). We define 




(2.1) 



Uijin) = Ui(n)Uj(n + i)Uj(n + j)U]{n) 



(2.2) 



b 



2N 2 



(2.3) 
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and take the large N limit with b held fixed. 

All computations were done on a L x x L y x L z periodic lattice with L x < L y = L z . 
One gauge field update of the whole lattice P] is one Cabibo-Marinari heat-bath update of 
the whole lattice followed by one SU(N) over-relaxation update of the whole lattice. The 
code was run on two clusters, one with 48 nodes and another with 31 nodes. The nodes in 
the cluster were simply used to generate more statistics using a parallel random number 
generator and generating independent configurations with the same set of parameters on 
different nodes. 

2.2 Determination of the critical sizes 

Given the lattice coupling b and lattice sizes L x and L y , the dimensionless physical sizes 
are defined as 

l xy = lim L xy /b t3A . (2.4) 

b— >oo 

The tadpole improved coupling j|] , 6 ta( j is defined as 

&tad = be(b) = K 12N * ^ Tr\Uij(n) + 17+ (n)]) (2.5) 

2.3 Lattice bulk transition 

Since the computations in this paper use the Wilson gauge action on the lattice, it is 
necessary to address the unphysical transition which is the extension of the Gross- Witten 
transition |7| in QCD2. The order parameter for this transition is the plaquette operator. 
The third order transition analytically computed in QCD2 remains to be true based on 
a numerical investigation in QCD3 ^ and the critical point is b = 0.43. This lattice 
transition does not survive the continuum limit and we will work with b > 0.43 throughout 
this paper in order to describe continuum physics. 

2.4 An order parameter 

An order parameter suitable for studying the phase transitions we are interested in is H 



Px,y,z — (Px,y,z) 

2L x LyL z 



1 2 

jfTrV x ,y, z (n) 



V x , y , z (n) = Ui(n + ml). (2.6) 



m=l 



The quantity P x ,y,z takes values in the range [0, 0.5] on any gauge field background and we 
choose the x, y and z directions on each configuration such that P x < P y < P z . 

Although this observable needs to be renormalized, we found it sufficient to work with 
the unrenormalized operator and we also did not have to smear the link variables. The 
eigenvalues, e l ° k ; k = 1, • • ■ , N, of the Polyakov loop operator, V x , y;Z (n), are gauge invariant. 
Px,y,z = 0.5 implies a uniform distribution of the eigenvalues of V Xy y )Z {n). A departure 
from P x ,y, z = 0.5 implies the presence of a peak in the distribution of the eigenvalues of 
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T- > x,y,z(iT') and a breaking of Zjy symmetry associated with this operator. There is no gap 
in the distribution of the eigenvalues of V x ,y,z when Zn is broken. 

3. Transitions in Polyakov loops 

All computations were done using N = 47. Having picked a lattice size L x x L y x L y , 
each run listed in Table |l| was a closed loop in the lattice coupling b. The third column 
in Table [l] shows the range of b and the step size in b. The fourth column, iVj, shows the 
number of thermalization sweeps at the two end points. Only one measurement was done 
per node at each b and the fifth column, iV&, shows the number of sweeps done at each 
intermediate b before the measurement. For example, the run on 3 3 lattice, started at a 
b = 0.5 and went up to a b = 2.5. A total of 2000 sweeps were performed at b = 0.5 and 
b = 2.5 and a total of 1000 sweeps were performed for all b in between 0.5 and 2.5. The 
step size in b was 0.05 and this code was run on the 31 node cluster with one measurement 
at each b per node. All values of b between 0.5 and 2.5 had two sets of measurements; one 
on the way up in b and one the way down in b. 
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Table 1: The parameters of all the runs used to study the transitions in Polyakov loops along 
with the results for the critical sizes. 



3.1 Details of the data analysis 

The plaquette as defined in ( |2.5[ ) was measured on all configurations and this was used to 
obtain the tadpole improved coupling, 6 t ad- Figure || shows the results for all three Polyakov 
loop observables as a function of g— for the data obtained on the 4 3 lattice (second row 
in Table |l|). The hysteresis is clear in both P x and P y for the lc-2c transition and it is 
seen in all the P x ,y,z for the 2c-3c transition. The critical size for the various transitions 
along with the error is obtained by locating the two points (one for upward direction and 
the second for the downward direction) where the error is largest in the observable that 
is broken. The vertical lines in Figure |l| shows the critical sizes along with the errors and 
these results are shown in the last three columns of Table |l|. 
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Within the 2c phase one sees a difference in P x and P y . But this is just a consequence 
of our choice of observable. Note that we have picked P x < P y on every configuration. If 
we assume two independent Gaussian random variables, a and /3, that have the same mean 
and variance, then one can show that the variables P x and P y defined as the minimum and 
maximum of a and (3 will be distributed such that 




Our data within the 2c phase is consistent with the above equation. 

We did not choose a range in b such that all transitions are seen on all L x , L y pairs 
since some of them were used only to investigate the lc-2c transition. But, we always 
picked a range such that the end points are in one of the four phases. 

3.2 Oc-lc transition 

Let us focus on the seventh column in Table |l| to study the confinement-deconfinement 
transition. The results on 3 3 , 4 3 , 5 3 and 6 3 show that the Oc-lc transition is physical since 
the critical size, f^S is the same on all four lattices within errors. The results here are 

d n, 

consistent with the older results presented in JlJ. We also studied the Oc-lc transition on 
4 x 5 2 and 3 x 4 2 and found that the critical size is independent of L y as expected. Figure || 
shows that the six results for the Oc-lc transition do scale properly and we estimate the 
continuum critical size to be l\ = 5.90(47). If we take the central value for the dimensionless 
string tension from ||, namely, yfa = 0.1975, then we get 

and this is consistent || with saying that ^ is the deconfinement temperature. 

3.3 lc-2c transition on L 3 lattices 

The physical size associated with the lc-2c transition, li- is expected to depend on l x , 
the temperature in the deconfined phase. We first estimate the critical size on lattices 
with L x = L y . We use the data on 3 3 , 4 3 , 5 3 and 6 3 . The four results show continuum 
scaling as can be seen from Figure [| and we conclude that the lc-2c transition exists in 
the continuum limit. We estimate 12(h) = 3.85(43). As mentioned before, lc phase is 
the deconfined phase. The system is a small finite box at zero temperature in the 2c 
phase. This transition occurs on a I 3 torus when the temperature is 1.53(21) times the 
deconfinement temperature. 

3.4 2c-3c transition on L 3 lattices 

We also investigated the 2c-3c transition on 3 3 , 4 3 , 5 3 and 6 3 . Here again, the four results 
show continuum scaling as can be seen from Figure ||| and we conclude that the 2c-3c 
transition also exists in the continuum limit. The transition size will depend on l x and l y 
when both of them are smaller than l z . But, we only estimate ^3(^3, 13) = 2.14(26) here. 
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Large N QCD on a very small torus, Z 3 , for I < 1% feels the size of the box and the 
temperature is high. Large N QCD is in a small box of size I at zero temperature if I > I3 
and it undergoes a phase transition into the deconfined phase when the box size is 1.80(30) 
times I3. 

3.5 Phase diagram for l x < l y < l z 

The single scale in the lc phase is l x which can also be thought of as inverse temperature in 
the deconfined phase. The 2c phase has two scales, namely the size of the two dimensional 
box l x and l y with l x < l y . If l y > £2(^)1 then the theory does not depend on l y and we 



are in the deconfined phase. We considered the special case of l x = l y in section 3.3. We 
extended this to the case when l x < l y . For this purpose, we considered the lattices listed 
in last six rows of Table |l[ The phase transition in P y is shown in Figure [|. There is an 
obvious dependence of the critical size I2 on l x . 

Figure |6| summarizes the various phases by focusing on the (l x , l y ) plane at l z =l y . The 
dependence of fo(lx) is shown using the shaded square points in Figure ^| The dashed line 
is a quadratic fit to the seven points and we note that ^(0) > 0. In order to get an overall 
picture, we have also shown the Oc-lc transition in Figure ^. The dotted line indicates that 
the Oc-lc transition does not depend on L for l y > l x . Figure |6| also shows the 0c, lc and 
2c phases for l x < l y < l z . For completeness, we have also shown the 2c-3c transition as 
seen on this specific (l x ,l y ) plane restricted to l x = l y . Like the lc-2c transition, the 2c-3c 
transition will also show a dependence on l x for l x < l y and the 2c phase will not reach 
the l x = l y line for l x = l y < fo- Furthermore, the 2c-3c transition curve will change as 
one changes the l z that defines the (l x ,ly) plane. We have not investigated these details 
pertaining to the 2c-3c transition in this paper. But, we should remark that the rest of the 
phase diagram does not depend on l z for l z > l y > l x . 



4. Conclusions 

Large N QCD in three dimensions on a / 3 continuum torus exists in four different phases. 
The theory is in the confined phase (0c) for I > 5.90(47) = l\ and physics does not depend 
on the box size. This critical size is the inverse of the deconfinement temperature, T c = ^, 
and the theory is in the deconfined phase (lc) for 1 < < 1.53(21). 

The system is in a finite box and feels the effect of temperature (3c phase) when 
/ < 2.14(26) = l 3 . The temperature has no effect if 1 < ^ < 1.80(30) (2c phase). The 
system goes into the deconfined phase if > 1.80(30). 

All phase transitions are most likely first order in nature. We have provided arguments 



for this scenario when going from lc-2c and 2c-3c. Previous results ||, 10 [ indicate that 
the deconfinement phase transition is also first order. 

The lc to 2c transition on a l x x l y x l z torus with l x < l y < l z depends on l x . The 
critical line is given by l 2 {l x ) = 0.56 + 1.08/^ - 0.059/^ and this is valid for < l x < 3.85(23) 
and it is independent of l z . This transition has been analyzed for one point on the (l x , l y ) 
plane in H and we are in agreement with the result in that paper. 
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Figure 1: Plot of P x ,y,z for the data in the second row of Table [T] showing all three transitions. 
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Figure 2: Plot of P x showing the Oc-lc transition. 
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Figure 3: Plot of P y showing the lc-2c transition on lattices with L x = L y 
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Figure 4: Plot of P y showing the 2c-3c transition on lattices with L x = L y . 
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